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Abstract 

Self-consistent calculations within the relativistic mean field the¬ 
ory (RMFT) were performed for 150 spherical even-even nuclei. The 
macroscopic part of the binding energy was evaluated by subtracting 
the Strutinsky shell corrections from the RMFT energy. The parame¬ 
ters of a liquid-drop (LD), like mass formula which approximates the 
RMFT results, were determined. The mass and isospin dependence of 
the RMFT mean-square radii constant for the neutron, proton, charge 
and total density distributions were estimated. The RMFT liquid- 
drop parameters and the radii constants are compared with similar 
results obtained with the Hartree-Fock-Bogoliubov calculations with 
the Gogny force (HFB-I-Gogny) and phenomenological models. 

RAGS numbers: 24.75.-G, 25.85.-w, 25.60.Pj, 25.70-z 


1 Introduction 

Self-consistent Hartree-Fock (HF) calculations with effective nucleon-nucleon 
forces of the Gogny and Skyrme type, or within the relativistic mean 
held theory (RMFT) Q are nowadays able to describe many features of 
nuclei. The theoretical results agree with the experimental data and the 
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masses, charge and neutron radii, electric multipole moments or energies 
of the lowest excited states are well reproduced even for nuclei beyond the 
stability line. It is probable that the presently accessible nuclei with the 
large neutron excess demand some revision of the parameters used in the 
traditional models, which have been adjusted to the smaller amount of data 
around the f3 stable nuclei. 

It is also interesting to compare the self-consistent prescriptions with 
other simpler models and see how they work for the nuclei close to the proton 
or neutron drip lines. The macroscopic-microscopic method with the liquid- 
drop (or droplet) model, using the Strutinsky shell correction and various 
kinds of the single-particle average potentials is of special interest because 
of its simplicity. Is it possible to extract the shell effects from the self- 
consistent energy and obtain an estimate of the macroscopic energy hidden 
in these models? This was already done successfully for the Skyrme [4] 
and Gogny [5] forces and now we would like to apply a similar method for 
the relativistic mean held (RMF). Nevertheless, one has to remember that 
the weak binding effects at drip lines are of purely quantal origin and the 
application of macroscopic-microscopic method could be questionable there. 

In Section II a short overview of the RMFT equations and parameters is 
given and the prescription for shell correction [6] is recalled. Moreover the 
liquid-drop formulae for the macroscopic energy and the root mean square 
(r.m.s.) radii of the proton, neutron and charge distributions are mentioned. 
In Section III the macroscopic part of the RMFT energy is approximated by 
a liquid-drop like formula. The liquid-drop parameters corresponding to the 
macroscopic part of the RMFT binding energy are compared with those of 
other theoretical and phenomenological models. In Section IV the RMFT 
root mean-square radii of 150 even-even spherical nuclei are approximated 
by the isospin-dependent formulae and compared with other experimental 
and theoretical estimates. The ratio of proton to neutron radius is of special 
interest due to the lack of experimental data for neutron radii. One can use 
this ratio to predict the neutron r.m.s. radius of a nucleus when its charge 
radius is known. At the end of the paper the conclusions are drawn and 
further investigations proposed. 
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2 Theory 

The single-particle level scheme obtained within the self-consistent RMFT 
calculation is used to evaluate the shell correction (-Esheii) to the binding 
energy: 

Esheii = X! ~ E , (1) 

occ 

where the sum runs over all occupied levels. All the single-particle levels 
up to the cut-off energy lying 15 MeV above the Fermi surface, are used to 
obtain the smoothed energy from the Strutinsky integral. We haven’t apply 
the newer prescription for the shell correction proposed in Ref. 0] to avoid 
the single-particle continuum effect because we wished to compare our results 
with those of reference obtained with the classical Strutinsky prescription. 
Nevertheless, the single-particle levels scheme for each nucleus, especially for 
the neutron-rich ones was carefully checked in order to take into account 
the proper number of single-particle states around Fermi surface. It was not 
our aim to estimate the position of the drip lines, but to obtain the average 
dependence of binding energy on A, Z number. For the detailed calculation 
of the binding energies of the nuclei close to the proton or neutron drip lines 
the use of the prescription of Ref. |0, H would be necessary. The Strutinsky 
smooth energy E is equal to 

E = 2 f ep{e)de . (2) 

J — OO 

The average levels density p{e) was obtained by the smoothing of the single¬ 
particle levels density p(e) = — e^) with the Gauss function multiplied 

by the 6 th order correction polynomial / 



The width parameter of the Gauss function 7 = 1.2hci;, with hoo = 
MeV , corresponds to the average position of the Strutinsky plateau 
in the shell corrections for the chosen sample of 150 spherical even-even nu¬ 
clei. The average single-particle levels density obtained with the RMFT is 
close to the results obtained in Ref. for the Gogny force. 

The macroscopic part of the binding energy is equal to the difference be¬ 
tween the self-consistently calculated RMFT energy (E^mpt) without pairing 
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interaction and the total (neutron and proton) shell correction 

rTiRMPT _ 771 rnn rnp (/i\ 

-^macr “ -^^RMFT - ^shell “ -^shell ' 

These quantities, evaluated for several nuclei with mass numbers A and 
isospins I = {N — Z)jA are approximated by the liquid-drop formula of 
Myers-Swiqtecki type 0 

-E'macr = ”&vol(l ” ^vol I^)A + 6surf(l “ «^surf l'^)A'^^^ + bcou\Z^A — C^Z"^ jA 

( 5 ) 

where 6coui is connected with the charge radius parameter by 6coui = 

3 2 /„ch 
5 ® /^O ■ 

The nucleon densities Pn(h') and Pp{r) obtained within the RMFT -|- BCS 
model could be used to evaluate the mean-square radii of the neutron or 
proton distributions. Here the question of the validity of that approach for 
the nuclei near drip line rises again, but the influence of pairing forces on the 
nuclear radius calculated in various models ( HFB -|- Gogny, LD -|- Woods 
Saxon ) is similar and does not interfere with the isotopic shifts, which have 
been measured 

{r^)q = J / J Pq{r)dV , q = {n,p}. (6) 


Knowing the mean square radius (r^) one can define an equivalent spherical 
sharp radius R using the following relation (r^) = , which arises directly 

from Eq. (j^) for the uniform density distribution pq = with Mq = 

{N,Z} and the volume conservation condition. In a rough estimate one 
usually assumes that R = tqA^A ^ takes the radius constant tq = 1.2 fm. 
However this formula turns out to be too approximate and it was proved 
in Ref. |]T| that a similar formula but using an isospin dependent radius 
constant, described the experimental data in a more satisfactory way. We 
have shown that the measured or calculated mean-square radii within the 
RMFT ||TI[ or HFB-|-Gogny p models could be accurately reproduced when 


the radius constant has the following form 


rQ = rooilR al + k/A) , (7) 

where roo, a, and n, are free adjustable parameters. The ratio of the proton 
to neutron root mean-square radii could be described by a formula similar to 
the one given above and could be used to predict the radius of the neutron 
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distribution when the charge radius is measured O, ^ . One has to note that 


this ratio does not depend on deformation in a first approximation since the 
density distributions of neutrons and protons are close to each other also for 
deformed nuclei. 


3 Binding energies 


The RMFT calculations with the NL3 set of parameters were performed for 
150 even-even nuclei between the proton and neutron drip lines which have, 
according to Ref. |^] a quadrupole moment almost equal to zero. They 
are: i 62 - 220 pp isotopes, = 50 with 

A E (86, 92), N = 82 with A E (122,164), and N = 126 with A E (174, 224) 
isotones and 30 other spherical nuclei along the (3 stability line. This choice of 
nuclei had already been used to estimate the shell effects by the HF method 
with the Gogny force . This set of representative spherical nuclei is larger 
than the sample of 30 deformed nuclei taken for the radii calculation within 
the RMFT in Ref. 
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We have used Nq = 20 shells and the oscillator length constant b = 
2.4 MeV of the harmonic oscillator as the basis when solving the self-consistent 
RMFT equations for fermions. At hrst the calculations were performed with¬ 
out taking into account the pairing residual interaction in order to evaluate 
the Strutinsky shell corrections, and then the experimental proton and neu¬ 
tron pairing energy gaps Ap, A„ were used to evaluate the r.m.s. radii and 
the potential energies in the RMFT -|- BCS model. This simplihed way of 
pairing correlation inclusion does not influence the values of radii significantly 
even for the nuclei near the drip lines. 


3.1 Liquid-drop parameters 

Fig. 1 shows the RMFT (solid lines) shell corrections in comparison with 
the results of Ref. [5] (dashed lines) obtained for the Gogny force. In the 
hrst panel of the multiplot one can observe the dependence on A of the total 
shell correction six groups of Ga-Th isotopes, in the middle i^gheii 

for three groups of N = 50, 82,126 isotones and in the r.h.s. panel for (3 
stable isotopes. 

The shell corrections obtained in both theoretical models are similar. 
They exhibit minima for the same magic numbers of one kind of nucleons and 
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differ from each other by no more than a few MeV. The RMFT estimates of 
the macroscopic part of the binding energy obtained by subtracting the total 
shell correction from the self-consistent RMFT energy (Eq. ^ for the above 
set of nuclei were htted by the liquid-drop formula (H), and the following set 
of parameters was obtained 


^RMFT 

= -15.19(1-1.66J2)R +16.81(1-1.21/2)yl2/3 +0.68-^-1.3— . 
MeV \ \ > ^ 1/3 ^ 

( 8 ) 

The r.m.s. deviation of the £t was equal to 1.97 MeV. In the Table I 
these RMFT estimates of the LD parameters are compared with the tra¬ 
ditional (MS-1967) Myers-Swi§,tecki liquid-drop formula 0. What is more 
Table I shows the modern phenomenological set (MS-2002) htted to 

0 


when using the microscopic 
In the last 
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presently available experimental masses 
(shell-|-pairing-|-deformation) energy corrections from Ref. 
column of Table I are given the results obtained in within the Hartree- 
Fock calculation with the Gogny DIS force |]I| which turned ut to be similar 
to these of RMFT. 

During the last 35 years, as seen in Table I, the liquid-drop parame¬ 
ters reproducing the experimental data have not changed very much. The 
macroscopic part of the binding energies obtained with the Gogny force 0 is 
described by the set of the LD parameters which approximates to the newest 
£t (MS-2002) of the LD parameters adjusted to the presently known exper- 

ill. 


imental masses 


The results obtained within the RMFT give smaller 


values of the volume and surface energies, while the charge radius constant 
corresponding to the Goulomb energy is equal to 1.264 fm and is substantially 
larger than its present phenomenological value (1.191 fm). By contrast, the 
RMFT estimate of the G 4 parameter, which is responsible for the charge dif¬ 
fuseness effect, is much closer to its phenomenological value compared with 
the Gogny’s one. The volume and surface dependence on isospin is weaker 
in the RMFT than the experimental one. The Gogny force gives a slightly 
stronger dependence of both energies than the phenomenological (MS-2002) 
one. 

We can compare the three models in Fig. 2. The results of MS-2002 
liquid-drop model and Gogny are subtracted from the macroscopic energies 
of the RMFT and shown for all the groups of nuclei in dependence on A. 
Since the RMFT macroscopic energy (solid lines) is the smallest, it gives the 
largest binding. The Gogny results are closer to the RMFT ones than to the 
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phenomenological (MS-2002) binding energy. 

The differences between the binding energies obtained with these three 
models reach even -30 MeV for isotope and isotone chains while for (3 stable 
isotopes they stay within -20 MeV. This is understandable because the NL3 
set of parameters of the RMFT was htted for nuclei close to the (3 stability 
line. The binding energies obtained with the Gogny force are closer to the 
liquid-drop estimatesthan these of RMFT. The isospin dependence of the 
binding energies is not well reproduced by either of the both models. 


4 Mean-square radii 


R is a known fact that the pairing correlations influence the density dis¬ 
tribution in nuclei. Therefore in order to evaluate the neutron and proton 
mean-square radii within the RMFT, we have to include the pairing forces. 
This was done in a simplihed way by inserting into the BCS equations the 
experimental proton and neutron energy gap between the ground state and 
the hrst excited two-quasiparticle state of even-even nuclei. The pairing en¬ 
ergy gaps are extracted from the experimental binding energies 
help of a three parameter formula proposed in Ref. 0 
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A, = - 1) - 2B{M,) + B{M, + 1)] , q = {n,p} , (9) 

where = (—1)-^'^ and Mq denotes nucleon number N for neutrons, Z for 
protons. When the BCS equations are solved the pairing correlations are 
added to the self-consistent mean held. 

The resulting r.m.s. radii for neutron and charge distributions as well as 
the ratio of proton to neutron radii are plotted in Figs. 3-5 for the three 
groups of isotopes, isotones and /3 stable nuclei. The RMFT radii can be 
easily reproduced by the isospin dependent formula (^, which corresponds 
to the sharp density distribution. The RMFT radius constants htted for the 
150 spherical nuclei are: 
for neutrons 


r” = 1.17(1 + 0.27/ + 3.38/71) fm , 

(10) 

for protons 


rl = 1.22(1 - 0.15/ + 1.51/71) fm , 

(11) 

and for charge distribution 


= 1.23(1 - 0.15/ + 2.47/71) fm . 

(12) 
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The r.m.s. deviation of each fit was smaller than 0.01 fm. 
([T0| - |T^ are very close to those obtained in Ref. 
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The estimates 
for the smaller sample 


(30) of deformed nuclei. This means that the deformation dependent function 


renormalizing the distributions to the sphere was properly chosen in Ref. [H 
and that the formulae ([l0|-|l^ adequately describe the radii constants , not 


only for the spherical but also for the deformed nuclei. 
This is also the case with the proton to neutron ratio 

^ = 1.04(1 - 0.38/ - 1.52M), 
r„. 


(13) 


which can be used to estimate the neutron radii with the help of the measured 
charge radius 

\Jrch - 0.64fm2 

“ 1.04(1 - 0.38/- 1.52/R) ^ ^ 

producing good agreement (with a slight tendency to overestimate) with the 
14 experimentally known neutron radii of Ref. |^. In contrast, a similar 
ratio obtained with the Gogny force in ||^ gives a slightly smaller neutron 
radius. Both groups of neutron radii for the 14 experimentally known data 
can be seen in Fig. 3 in dependence on the reduced isospin I. 

In Fig. 4 the differences between the charge radii predicted by the RMFT 
and the Gogny model (solid lines) are compared with the experimental data 
from which also the Gogny radii are subtracted (crosses). One can 


see that the agreement of the RMFT results for the charge radii with the 
experimental data is even slightly better than that of Ref. obtained with 
the Gogny force. 

In Fig. 5 the proton to neutron radius ratio obtained in the RMFT 
(solid lines) and with the Gogny force (dashed lines) is compared with the 
experimental data (crosses) |^, ^ for the three groups of isotopes, isotones 
and f3 stable nuclei. 


The parameters of formulae ([T0| - |l2D obtained for various theoretical mod¬ 
els are compared in Table II with the ones htted to the experimental data 
T9| and in Ref. 
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10] for charge radii. 


Both self-consistent theoretical models give similar estimates of the roo 
parameter of neutron, proton, and charge radii . The isospin dependence 
of the r.m.s. radii is slightly different in the two models. The k/A term, 
important for the light nuclei, shows some differences as well. 














































5 Conclusions 

The following conclusions can be drawn from our investigation: 

1. The shell corrections obtained in the RMFT with the NL3 set of pa¬ 
rameters and within the Hartree-Fock mean held calculation with the 
Gogny DIS force are similar. 

2. The volume and surface parts of the binding energy in the RMFT 
are smaller than the corresponding energies obtained with the Gogny 
model [Q as well as than those of the liquid-drop model htted to the 
experimental masses P, |^. 

3. The isospin dependence of the volume and surface term obtained within 
the RMFT is too small in comparison with the phenomenological liquid- 
drop model. 

4. The mean-square radii of the proton, neutron and charge distributions 
are similar in the Gogny and RMFT models. 

5. The RMFT ratio of proton to neutron radii, used to predict the neu¬ 
tron radii when the charge radius is known, gives the estimates within 
experimental error bars for all 14 experimentally known neutron radii. 

Similar effects for deformed nuclei with various sets of RMFT parameters 
will be investigated soon. 
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Table 1: The macroscopic energy parameters. 


parameter 

unit 

MS-1967 

MS-2002 

RMFT 

Gogny 

^vol 

MeV 

15.667 

15.848 

15.185 

15.649 

^vol 

- 

1.790 

1.848 

1.657 

1.916 

&surf 

MeV 

18.560 

19.386 

16.811 

18.928 

^surf 

- 

1.790 

1.983 

1.209 

2.108 

.«ch 

' 0 

fm 

1.205 

1.190 

1.264 

1.188 

C'4 

MeV 

1.211 

1.200 

1.299 

2.015 
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Table 2: The radii parameters. 



Then. 

Gogny (150 sph.n.) 

RMFT (30 def.n) 

RMFT (150 sph.n.) 

neutrons 

m 

i 

El 


roo 

1.17 

1.17 

1.17 

1.17 

a 

0.16 

0.12 

0.25 

0.27 

K 

3.85 

3.29 

2.81 

3.38 

protons 

roo 

18 

1.22 

1.21 

1.24 

1.22 

a 

-0.17 

-0.14 

-0.16 

-0.15 

K 

1.78 

1.83 

0.65 

1.51 

charge 

roo 

0(0) 

1.24(1.25) 

1.22 

1.24 

1.23 

a 

-0.19(-0.25) 

-0.15 

-0.15 

-0.15 

K 

1.65(2.06) 

2.32 

0.58 

2.47 

ratio 

roo 

00 

1.03 

1.04 

1.05 

1.04 

a 

-0.36 

-0.27 

-0.36 

-0.38 

K 

-1.33 

-1.12 

-3.15 

-1.52 
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Figure 1: The total shell corrections obtained within the RMFT (solid lines) 
and with the Gogny force (dashed lines) in dependence on the mass number 
A. The three parts of multiplot show the shell corrections of Ca-Th isotopes, 
for the N = 50, 82,126 isotones and for the (3 stable nuclei respectively. 
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Figure 2: The comparison of macroscopic energies obtained within the RMFT 
and with the Gogny ^ force and the liquid-drop energies []^ in dependence 
on A. The differences between the RMFT macroscopic parts of the binding 
energies (solid lines) and the phenomenological (MS-2002) (solid lines) 
estimates are compared with the corresponding differences between RMFT 
and the Gogny macroscopic energies (dashed lines) for the isotopes (l.h.s.) 
isotones (middle) and /3 stable nuclei(r.h.s.). 
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Figure 3: The 14 experimentally known neutron radii from Ref. (crosses) 
are compared with the RMFT predictions (spheres) and the estimates done 
with the Gogny model (stars) in dependence on the reduced isospin / = 
(N-Z/A). 
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Figure 4: The charge radii predicted by the RMFT related to the results 
of the Gogny model (solid lines) are compared with the experimental data 
||T8| from which the same reference of Gogny charge radii is also removed 
(crosses). The three panels correspond to the isotopes (up) isotones (middle) 
and f3 stable nuclei (down). 
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Figure 5: The proton to neutron radius ratio obtained in the RMFT (solid 
line) and with the Gogny force (dashed lines) is compared with the similar 
ratio obtained on the basis of the experimental data (crosses) [|T^ |T^. The 


three parts of multiplot show the calculated root mean-square radii ratios 
for the three groups of isotopes (l.h.s.), isotones (middle) and f3 stable nuclei 
(r.h.s.) §. 
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